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ABSTRACT 

The Dual Imaging Cerenkov Experiment (DICE) was designed and operated for mak- 
ing elemental composition measurements of cosmic rays near the knee of the spectrum 
at several PeV. Here we present the first results using this experiment from the mea- 
surement of the average location of the depth of shower maximum, < X max >, in the 
atmosphere as a function of particle energy. The value of < X max > near the instrument 
threshold of ~ O.lPeV is consistent with expectations from previous direct measure- 
ments. At higher energies there is little change in composition up to ~ 5PeV. Above 
this energy < X max > is deeper than expected for a constant elemental composition 
implying the overall elemental composition is becoming lighter above the knee region. 
These results disagree with the idea that cosmic rays should become on average heavier 
above the knee. Instead they suggest a transition to a qualitatively different population 
of particles above 5 PeV. 
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1. Introduction 

Observations of cosmic rays have presented several 
puzzles which have been resolved by detailed mea- 
surements of composition, both elemental and iso- 
topic. The present picture of the production of the 
bulk of cosmic rays by diffusive shock acceleration in 
supernova remnants is attractive since it appears to 
be the only energy source powerful enough to sus- 
tain this particle population in our galaxy. There is 
also a strong theoretical link between the observed 
energy spectra of the sources and the predictions of 



strong shock acceleration, ( 


Axford et al. 1977 




Bell 


1978. 


Blandford and Ostriker 1978, 


Krymsky 1977). 



However, there is still no direct evidence that the 
bulk of cosmic rays are produced by this mecha- 
nism. Searches for the TeV gamma emission which 
should be produced by the high energy nuclei at these 



sites are in progress, (see e.g. Lessard et al. 1995). 
The supernova remnant shock acceleration in its most 
straightforward form produces a maximum energy for 
particles accelerated by this mechanism to be close 



to ZxO.lPeV flLagage and Cesarsky 198% This limit 
is determined by the average acceleration rate and 
the mean lifetime of the strong shock in a typical 
supernova event. The 'knee' in the cosmic ray flux, 
where the energy spectrum becomes steeper around 
several PeV, could well be re lated to this theoretical 
limit. Some workers (see e.g. Axford 1994 ) have sug- 
gested that cosmic rays above the knee are produced 
by the secondary acceleration of high energy parti- 
cles from supernova by weaker shocks in the Galaxy. 
This would result in an increase in the average mass 
of cosmic rays through the knee. If the cosmic rays 
from shock acceleration in Galactic supernovae are 
more strictly limited to the knee energy then a sec- 
ond population of particles must be present to pro- 
duce the fluxes of cosmic rays observed at higher en- 
ergies, possibly of extragalactic origin as suggested by 
Prothcroc and Szabo 1992. As with other properties 



of cosmic rays a measurement of composition in this 
energy range can be used to detect these changes and 
provide the information needed to investigate the ori- 
gins of cosmic rays above the knee. 

While direct detection above the atmosphere re- 
mains the most desirable method for composition 
determination, the fluxes of particles in this energy 
range are so small that a space instrument large 
enough to collect a good statistical sample up to 
lOPeV has not yet been flown. Indirect detection 



of cosmic rays, through measurements of air showers 
produced in the atmosphere, can easily supply suf- 
ficient collecting power but these methods generally 
have substantially reduced mass and energy resolu- 
tion. The properties of the originating cosmic ray 
nucleus becomes clouded by the huge number of in- 
teractions in the air shower. The reconstruction of the 
mass of the incoming nucleus from measurements of 
the shower distributions at ground level is also subject 
to systematic errors resulting from imperfections of 
the Monte Carlo modelling of air showers. DICE is a 
ground based air shower experiment which is designed 
to have as little reliance as possible on the details of 
the Monte Carlo simulations and to have the capabil- 
ity of comparison of results with direct measurements 
at O.lPeV to provide an assessment of overall system- 
atic errors. 

2. Experiment Description 



The two DICE telescopes (see Boothby et al. 1995 
Boothby et al. 199^ ) are located at the CASA- 



MIA site in Dugway, Utah (described in Borionc 
ct al. 1994). They each consist of a 2m diameter 



f/1.16 spherical mirror with a focal plane consisting 
of 256 close packed 40mm hexagonal photomultiplicrs 
(PMTs) to provide ~ 1° pixels with a field of view of 
16° x 13.5° centered about the vertical. The two tele- 
scopes are on fixed mounts separated by 100m located 
near the center of the CASA site. 

Cosmic ray events within the field of view produce 
a focal plane image at the photomultiplicr cluster 
which is an intensity pattern of Cerenkov light com- 
ing from the air shower. When the direction of the air 
shower and the distance of the shower core from the 
telescopes are known, simple geometry can be used 
to reconstruct the amount of light received from each 
altitude of the shower. The amount of Cerenkov light 
produced at each altitude is strongly correlated with 
the number of electrons at that atmospheric depth. 
This is used to estimate the electron size as a func- 
tion of depth in the atmosphere from which X max can 
be determined. This procedure is essentially geomet- 
rical and is independent of Monte Carlo simulations 
except for calculations which determine the angular 
distribution of Cerenkov light around the shower axis. 

Here we report on the analysis of shower maxima 
from the images in the DICE clusters using the inde- 
pendent shower geometry established by the CASA 
air shower array. The accuracy of the CASA geometry 
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for each shower in the knee region is < 5m for shower 
core location on the ground and < 0.6° for shower 
trajectory direction (Borione et al. 1994). The ab- 
solute alignment of the DICE clusters relative to the 
vertical has been determined to an accuracy of ~ 0.1° 
by detection of transits of stars across the DICE aper- 
tures. 

The events which trigger both DICE telescopes and 
the CASA array have data which are recorded sepa- 
rately and later matched by comparing the times of 
each event. The trigger rates are typically 0.1 Hz and 
the time information is accurate to ~ 1ms, the frac- 
tion of mismatched events in the raw data is <1%. 
During analysis these mismatches are further sup- 
pressed by (i) the requirement that shower images 
in DICE lie in the direction of the CASA event, (ii) 
the reconstructed shower maxima from the two DICE 
telescopes using the CASA geometry are in agree- 
ment. These analysis cuts further suppress the mis- 
matches by a factor of ^1000, resulting in an overall 
background level of ~ 10~ 5 . 

3. Data Collection and Analysis 

DICE has been operational at the Dugway site 
since summer 1994, this analysis comes from the pe- 
riod July 1994 - December 1995. Data is collected 
during clear, moonless nights. The events analyzed 
here come from 914 hours of operation, correspond- 
ing to ~ 7% livetime during this period. The effec- 
tive geometrical acceptance aperture of DICE with 
the data cuts used is ~2000 m 2 sr making the over- 
all collecting power at high energies near 80,000 m 2 
sr days. This allows a statistically valuable measure- 
ment of composition up to energies of ~10 PeV. The 
trigger threshold and logic of the DICE PMT clusters 
produces some trigger inefficiency for showers below 
^800TeV, this introduces a small bias in this region 
which is discussed later. 

For each collected event the Cerenkov images are fit 
by translating the PMT amplitudes into a sequence 
of Cerenkov luminosities and hence shower electron 
sizes versus atmospheric depth (g/cm 2 ) by using the 
shower geometry from CASA and the distribution of 
Cerenkov light about the shower axis predicted from 



the CORSIKA air shower Monte Carlo (Knapp et al 
1996). These shower development curves are fit to 



a longitudinal development function from which the 
shower maximum, size and width are determined. In- 
dependent fits are made for each of the two DICE 



clusters so that a consistency check can made between 
the two sets of shower fit parameters. For the shower 
maxima the difference distribution between the fits 
from the two sites has an rms width of ~65g/cm 2 , 
after correcting for correlations associated with the 
uncertainties in the geometry the dual site rms reso- 
lution for a single event is ~45g/cm 2 . Showers which 
have cores located equidistant from each cluster pro- 
duce amplitudes of Cerenkov light which should be 
similar in each cluster. The difference distribution of 
these events shows an rms width of ~26% implying 
a dual site resolution of ~13% for the Cerenkov size. 
The energy of each event is estimated from the mean 
of the independent values derived from the light in- 
tensity at each cluster. The amount of Cerenkov light 
collected at each cluster is corrected by a factor of 
B x exp(r/ro) where B is a normalization factor, ro is 
a scale length, and r is the distance from the cluster to 
the core. The values of r and B are determined from 
CORSIKA. These do not have a large dependence on 
the mass of the primary nucleus for the core distances 
of 100m< r <225m used in this analysis. In this anal- 
ysis the value r =82.5m is used and the value of B is 
taken to have a single value for all events. This intro- 
duces a mass dependent error in the energy estimate 
which is ~10%. The largest uncertainty in the energy 
estimate is in the conversion from PMT pulse height 
to Cerenkov photon density since a number of fac- 
tors including mirror reflectivity and PMT quantum 
efficiency need to be accurately known. This overall 
systematic error in the energy scale is estimated to be 
±20%. 

The following data cuts are made: (i) The core of 
the shower lies at a distance 100m< r <225m from 
both DICE clusters. (ii)The fits of the longitudinal 
developments in each DICE telescope have reduced 
X 2 < 3. (iii) The X max from the two sites agree within 
75g/cm 2 . (iv) The shower direction is within 5° of the 
vertical. The measurements of X max and energy from 
the two clusters are averaged to give overall values for 
each shower. Values of X max are binned versus en- 
ergy using the boundaries shown in Table |l|. To gain 
some confidence for the accuracy of the energy scale 
and the overall detection efficiency estimates an en- 
ergy spectrum of cosmic rays has been derived from 
the data by calculating absolute fluxes. The dominant 
efficiency correction for this work is the trigger scheme 
for each photomultiplier cluster. A simulation which 
includes the trigger logic and air shower character- 
istics derived from samples of air showers produced 
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by CORSIKA is used to derive the various efficien- 
cies and the acceptance geometry of DICE. The ex- 
posure livetime can be directly derived from the data 
since at each minute of operation the cluster electron- 
ics records the livetime for that period. The cluster 
deadtime is small, typically <1%. A combination of 
the efficiencies, geometry and livetime is used to de- 
rive the fluxes of particles by a matrix deconvolution 
technique using the monte carlo. 

4. Results 

The absolute fluxes of cosmic rays obtained for 
these data are shown in Figure [I], note the fluxes 
are multiplied by an overall factor of E 2 75 . Also 
shown for comparison are results from some previ- 
ous air shower experiments which use electron sizes 
at ground level to estimate the energy. The DICE 
data energy spectrum has been calculated assuming 
a mixed composition model of the type discussed by 



3 wordy 1995. This model provides a composition at 
IPeV of -40% H, -30% He, and -30% heavier nuclei 
and at lOPeV -20% H, -20% He and -60% heavier 
nuclei. Also shown are energy spectra reconstructed 
under the assumption that either (a) the particles are 
all protons or (b) particles are all iron nuclei. As can 
be seen these do not change the fluxes or slopes sig- 
nificantly. There is a distinct change in slope near 
— 4PeV which is consistent with the existence of the 
'knee' in this region. Since the model used to recon- 
struct the energy spectrum is heavier at high energies 
than suggested by the values of < X max > measured 
by DICE, this could introduce a systematic bias which 
might enhance the derived fluxes above — 5PeV. This 
enhancement is at maximum 10%, less than the sta- 
tistical precision of the fluxes in this region. The pos- 
sible effect of a miscalibration of the energy scale is 
shown by the length of the vertical bar to the left 
of the picture. This represents the overall system- 
atic uncertainty in flux from the DICE measurements 
corresponding to an energy scale shift of ±20%. The 
size of the trigger inefficiency correction becomes large 
below — 400TeV making fluxes derived from events in 
this region more uncertain than this systematic limit. 
The overall DICE trigger efficiency is a combination 
of the cluster trigger efficiency and the total light in- 
tensity produced at the cluster. DICE is therefore 
capable of triggering with high efficiency on bright 
events with energies >2PeV out to the maximum ex- 
tent of the data cuts at r = 225m. The light from 
lower energy particles showers is not intense enough 



to always provide triggers at this distance, the trig- 
ger simulation is used to calculate the amount of data 
lost by this effect. At 500TeV the overall detection 
efficiency is —50%. 

The variation of < X max > with energy is given in 
Table [l]. The third column shows the raw values and 
the next column shows the values of the correction, 8, 
for trigger and geometry bias. The trigger bias arises 
from the fact that near threshold protons are more 
likely to trigger DICE than iron nuclei because they 
produce more Cerenkov light. The size of this cor- 
rection is derived from the trigger simulation which 
incorporates information from CORSIKA on the vari- 
ation of the intensity of light with incoming particle 
energy and mass. The geometry bias of this correction 
is associated with the finite pixel resolution of DICE. 
At high altitudes the shower structure becomes dif- 
ficult to determine because the shower development 
occurs on angular scales which are comparable to the 
DICE pixel resolution. The combination of these ef- 
fects produces an overall correction which is relatively 
large for the first energy bin but which becomes small 
at higher energies. These data are plotted in Figure |[ 
the vertical errors are statistical, the horizontal cor- 
respond to the systematic error in energy estimate 
discussed above. Also shown is the results of direct 
measurements at O.lPeV converted by the CORSIKA 
Monte Carlo to a value which corresponds to the ex- 
pected mean shower maximum for this composition. 
The small dotted lines show how the mean shower 
maximum should behave for pure protons or iron nu- 
clei based on CORSIKA. The DICE results lie be- 
tween these extremes and agree well with the extrap- 
olation from direct measurements at the lower ener- 
gies, shown as the large dashed line. There is little 
apparent change in mass composition through the en- 
ergy range of most of these data. The indications are 
that the composition becomes lighter above — 5PeV 
although the statistics are limited. The comparison 
of the measured distribution widths of X max with the 
widths of simulations for pure protons and pure iron 
nuclei including the DICE detector response are given 
in Table |[ These are consistent with the movement 
of the mean shower maximum with energy. 

The first results from DICE show no evidence for 
an increase in the average mass of cosmic rays in 
the knee region. These lower energy results of DICE 
are consistent with recent measurements which use 
the intensity distribution of Cerenkov light at ground 
level to show no change in composition up to IPeV 
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( Plaga et aL 199q ). An increase in mass has been 
expected somewhere near the knee because any pro- 
cess which involves a steepening of magnetic rigidity 
(momentum/charge) spectra of the particles, either in 
the sources or elsewhere, naturally leads to an over- 
all increase in the mass of a population of particles 
measured by energy/particle. It is difficult to escape 
the conclusion from Figure |^ and Table ^ that a sim- 
ple rigidity steepening of this type is probably not 
responsible for the 'knee' of the cosmic ray spectrum. 
Although these measurements indicate a transition to 
a more proton-rich source of particles, in the present 
results this is only a ~ 3cr effect which should be veri- 
fied by further measurements. If confirmed, this tran- 
sition suggests we are reaching an energy scale where 
conventional supernova shock acceleration and prop- 
agation through the galaxy is no longer the dominant 
source. We may be beginning to see a population of 
much older particles confined in the halo of our galaxy 
or possibly of extragalactic origin. 
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Table 1 

DICE MEASUREMENTS OF < X max > 



Energy bin 
(PeV) 


Events 


< X max > raw 
(g/cm 2 ) 


(5 a 
(g/cm 2 ) 


Median Energy b 
(PeV) 


(g/cm 2 ) 


0.2 - 0.4 


6132 


474 


-12 


0.25 


462±2 


0.4 - 0.7 


3708 


479 


-1 


0.5 


478±2 


0.7- 1.5 


2279 


501 


+2 


1.0 


503±2 


1.5 - 3.0 


740 


534 


+3 


2.0 


537±3 


3.0 - 5.0 


212 


549 


+3 


4.0 


552±4 


5.0 - 8.0 


70 


587 


+3 


6.0 


590±10 


8.0 - 15.0 


37 


628 


+3 


10.0 


631±14 



a This trigger correction is discussed in the text 

b The energy scale has an estimated systematic error of ±20% 



Table 2 

dice measurements of the widths of x max compared with simulations for protons and iron 

NUCLEI 



Median Energy 3, 
(PeV) 


DICE a(X max ) 
(g/cm 2 ) 


Sim. a (protons) 
(g/cm 2 ) 


Sim. a{iron) 
(g/cm 2 ) 


1.0 


73±1 


87 


66 


2.0 


66±2 


83 


49 


4.0 


58±4 


82 


41 


6.0 


82±12 


81 


40 


10.0 


85±20 


79 


40 



a Only energies where the trigger bias is small are included 



6 



REFERENCES 

Amcnomori, M. et al. 1996, ApJ, 461, 408 

Axford, W.I., Leer, E., Skadron, G., 1977, Proc. 15th Int. Cosmic Ray Conf. (Plovdiv) , 132, 11 

Axford, W.I., 1994, ApJS, 90, 937 

Bell, A.R. 1978, MNRAS, 182, 147 and 443 

Blandford, R.D. and Ostrikcr, J.P., 1978, ApJ, 221, L29 

Borionc, A. et al., 1994, Nucl. Instr. Meth., A346, 329 

Boothby, K., Chantell, M., Green, K.D., Kieda, D.B., Knapp, J., Larsen, C.G., Swordy, S.P., 1995 Proc. 24th Int. 
Cosmic Ray Conf. (Rome), 2, 697 

Boothby, K., Chantell, M., Green, K.D., Kieda, D.B., Knapp, J., Larsen, C.G., Swordy, S.P., 1997 Nuclear Physics 
B (Proc. Suppl.) 52B, 166 

Knapp, J., et al.(CORSIKA 4.50 with GHEISHA, VENUS), 1996 Forschungszcntrum Karlsruhe, Report FZKA 
5828 

Krymsky, G.F., 1977, Doklady Akad. Nauk. USSR, 234, 1306 
Lagage, P.O. and Cesarsky, C.J., 1983, A&A, 125, 249 
Lessard et al. 1995 Proc. 24th Int. Cosmic Ray Conf. (Rome), 2, 475 
Nagano, M. et al. 1984, J. Phys. G, 10, 1295 

Plaga, R. et al. 1995 Proc. 24th Int. Cosmic Ray Conf. (Rome), 2, 693 
Prothcroe, R. J. and Szabo A. P. 1992, Phys. Rev. Lctt.69, 2885 

Swordy, S.P., 1994 Proc. 23rd Int. Cosmic Ray Conf. Invited Papers, World Scientific: Singapore, 254 
Swordy, S.P., 1995 Proc. 24th Int. Cosmic Ray Conf. (Rome), 2, 697 



This 2-column preprint was prepared with the AAS IAT^X macros v4.0. 



7 



10° - 



c 
'o 

■4-* 

l_ 

Q) 

O 

c 

_o 

£ 

0) 

-M " 

OT 

>% 

CO 



□ 



□ □ 



□ □□□□□ ____ A 4____^ (b) 




J I I I I L_l_ 



_l I I I I I L_l_ 



10 



10 



Energy(PeV) 



The absolute fluxes of particles measured by DICE as filled circles compared with previous measurements 

1996| ). Note the fluxes are mu ltiplied by a f actor of 

Curve 



Fig. 1 

(triangles - Nagano et al 

(Energy) 2 . These fluxes were derived assuming the mixed elemental composition from Swordy 1995 



1984, squares -Amenomori et al. 



(a) shows the fluxes produced by assuming an all-proton composition, curve (b) shows fluxes assuming an all-iron 
composition. The vertical bar on the left represents the systematic variation of the flux levels because of uncertainty 
in the overall energy scale. The DICE fluxes have power law indices of -2.55±0.01 below 4PeV and -2.92±0.1 above 
4PeV 
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Fig. 2. — The variation of the mean depth of shower maximum (< X max >) with p article energy measured by 
DICE. The star is a conversion of an assembly of direct measurements ( [3 wordy 1994 ) to the expected < X max > 
using the CORSIKA Monte Carlo - the long dashed line shows how < X max > scales with energy for a constant 
elemental composition. The upper dotted line is the < X max > expected for all-proton composition and the lower 
dotted line is the value expected for all-iron composition. 
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